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Abstract. This paper describes some generalities about spectro-interferometry and 
the role it has played in the last decade for the better understanding of circumstellar 
matter. I provide a small history of the technique and its origins, and recall the basics of 
differential phase and its central role for the recent discoveries. I finally provide a small 
set of simple interpretations of differential phases for specific astrophysical cases, and 
intend to provide a "cookbook" for the other cases. 



1. Matter with movement which matters 

When asked to present "circumstellar matter studied with interferometry", I was put in 
a difficult situation, given the enor mous amount of literature on such a subject (see e.g. 
figure 2 in iLe Bouquin et alJl2009h . In the 80's (when I was born) this would have been 



much easier than today, as only a handful of stars had been observed with the very few 
interferometers existing. In-between, large collecting area interferometers and general- 
user instruments have appeared, making the former observing achievements almost as 
easy as getting a spectrum with a spectrograph. 

When speaking of matter around a central object, one can think of gas and/or dust, 
but one cannot separate its composition from its movement around that object. Indeed, 
a static sphere of dust would not wait very long until falling freely onto the centre 
of the system. Therefore, to make a stable system, one would need to impulse some 
movement to that sphere in order to counteract the gravity force: 

• One would need, for example, to impulse rotation to this sphere in order to pre- 
vent it from free fall. But in this case, angular momentum effects would flatten it 
to a disk-like structure. 

• If expansion (or contraction) would be impulsed to the sphere, then nothing 
would change the shape of the sphere, but the gravity force would slow down 
(or accelerate) that sphere. 

• More complex physical effects could also occur: 

- like in rotating and contracting viscous disks, where the loss of angular 
momentum would force the energy to escape from the system, either via 
heating, or via kinetic energy (materialized by jets). 

- such complex distribution of matter (disk+jet) would form a brake to an 
explosion of the central object (like in novae or supernovae), shaping the 
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remnant neb ula into bipolar outflows dWeigelt et all 120071 : IChesneau et al. 



120071120111) . 

- or like in turbulent flows, where the global structure of the flow (rota- 
tion/expansion) is respected, but not the local structure. 

One can see here that nothing about the mass or the nature of the central object was 
hypothesized, and therefore, studying the kinematics of matter around a central object 
applies to many classes of objects with the same techniques. These objects range from 
the accretion disk around the white dwarf of a nova system, up to the accretion disk 
around the supermassive black hole in an active galactic nucleus (when sufficiently far 
away to neglect relativistic effects). Such accretion disks can als o be found around a 



wide range of YSOs, f rom T Tauri stars (Ivan Boekel et al] 120041) up to very massive 
stars dKraus et alj|2010h . 

Studying such velocity fields has been made for decades with spectroscopy, but 
classical spectroscopy lacks the very important combination of velocity (spectroscopy) 
and position (astrometry) for all these compact objects. Interferometry provides this 
missing astrometric information, making possible that combination. Therefore, I will 
try to restrain in this paper to the study of the combination of the distribution of matter 
and its movement around a central object. 



2. The origin of spectro-interferometry 
2.1. At the beginning 

Speckle interferometry dLabeyridl 19701) revolutionized high-angular resolution imaging 
in astronomy. However, in the first times of this technique, no phase measurement was 
attainable, making the obtained high-angular resolution information very difficult to 
interpret. Indeed, with only an access to the amplitude of the Fourier components of 
the image of the object, this technique was only able to provide autocorrelations of 
the object's i mage. One w orkaround of this limitation, called speckle holography, was 
introduced by lGezari et al.l(ll97ll) . with a spectacular application to the r esolution of the 
centr al object in R136a, in the core of the young stellar cluster 30 Dor (Weigelt & Baier 



1985). However, it needed an unresolved reference star inside the anisoplanetism angle 

of the atmosphere in order to work. 

Until the advent of triple correlation techniques dHofmann & Weigeltlll993l) . sev 



eral ideas were proposed to recover the phases of the objects, by using the ir deter- 
ministic change as a func tion of spatial frequen cies dKnox & T hompsonl l 19741) . or as a 
function of wavelengths dKoechlin et al.|[l979l) . However, the first to propose to com- 
pute phase differences between adjacent wavelengths in speckle interferometry as an 
observable is iBeckersI dl983l) . with a direct potential application to close double stars 
research. In that work, the author proposed to use the properties of the cross-correlation 
of speckle patterns on the Multi-Mirror Telescope (MMT) to compute wavelength- 
dependent photocenter displacements of the telescope Point Spread Function (PSF). 
This displa c emen t is linked directly to the nature of the object, here binary stars. In 
lAime et all dl984l) . the idea of Beckers was first tested on-sky. At this time, the cross- 
correlation technique was used on single-dish telescopes (the MMT being considered 
here as single dish). 

Later, this technique was applied to long-baseline int erferometry on the Grand In- 
terferometre a deux Telescopes (GI2T), with some success dStee et al.ll995klVakili et all 
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1997111998b Berio et al. Il999h . Here, the cross-correlation technique was used to com- 
pute a phase as a functi on of wavelen gth on long-baseline interferometry data. 

Spectroas trometry [Baileyl (1 19981) was later invented, with the same purpose as the 
Beckersl (l983) idea, i.e. study the photo-center shift as a function of wavelength, but 



with a different algorithm (in the image space rather than Fourier-space). 



2.2. Today's spectro-interferometry 

Today, spectro-interferometry has bee n popularized by th e Very Larg e Telescope Inter- 
ferom eter (VLTI), with both AMBER dPetrov et al.ll2007l) and MIDI dGraser & Leinert 
1998) instruments, joined by the Navy (Prototype) Optical Interferometer (NPOI, now 
NOI, the " prototype" asp ect having disappeared n ow) in its diffe rential phase referenc- 
ing mode (ISchmitt et al.ll2009l) , VEGA/CHARA dMourard et"aDl2009h . and Keck used 
in self-referenced phase mode (Wo illez et al. II2012J). 

The results coming from these instruments are significant, among which the di- 
rect detections of the formation zones of se v eral ions, molecules, or dust types in 
young stellar objects ( van Boe kel et all 12004 iMalbet et all 120071 : iTatulli et all 12007 , 
20081 : lEisner"& Hiflenbr andTkoill). the first direct detection of a Keplerian-rotating disk 



aroun d a Be st ar (Meilland et al. 2007b) follow ed by many simi lar results on other Be 
stars (Meillan d rtali2007al : ICarciofi et al.ll2009l : IStefl et al.l2012l :lMeilland et al. 2012), 



the C O flow motion in the turbulent atmosphere of Betelgeuse (Ohnaka et al .l 1200' 



2008 



2011) , or the c haracteriz ation of binary system s (IMillour et all 120071 : Meillan d et al. 



Schmitt et al.ll2009l : iMeilland et al.ll2011 



More recently, a few attempts have been made to include the differential phases 
into an image-reconstructio n process, either ma king assumptions on the target's reso- 
lution by the interferometer ( Schmitt et al.ll2009l). or using radio-like image reconstruc- 
tion methods in a more general case (IMillour et al.ll201 lMOhnaka et al.ll201 1 ). These re- 
cent attempts illustrate the fact that differential phase do add more information about the 
observed object than the one contained in classical interferometric observables alone, 
namely squared visibilities and closure phases, and therefore its huge potential for im- 
proving the image reconstructions with optical interferometers. 



3. How to get it? 

The main idea behind computing interferometry observables is to get the invariant part 
of the wandering fringes of an interferometer, as seen through the atmosphere and a 
spectrograph. 

In the classical methods, this fringe wandering can be nulled in specific cases, 
where it has been demonstrated that no (or few) atmospheric effects affect the chosen 
observables. 

For example, one can get the Fourier amplitudes using speckle techniques (i.e. 
using the power spectrum instead of directly the Fourier transform of the fringes). This 
calculation provides V 2 (squared FT amplitude) as a function of wavelength. 

However, for computing the Fourier phases, things get more complicated: 



• One way is to use the closure relation when usin g three telescopes, the very same 
way as for bispectrum speckle interferometry dHofmann & Weigeltlll993 ). One 
phase (out of three) can therefore be extracted, independent from any instrumen- 
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tal or atmospheric effect. This phase is called "closure phase", but when using 3 
telescopes, it misses 2/3 of the available phase information. 

• Another way is to have at his disposal a model of the variable part of these 
wandering fringes, that can be subtracted from the raw Fourier Transforms before 
averaging. One does not need to have a model of the observed object, but rather 
a model of the atmosphere and the instrument to correct for these effects and 
calculate a phase (called "differential phase"). A by-product of this method is that 
one can get also the Fourier amplitude (called "linear visibility", or "differential 
visibility"). 

This last point is what is called differential interferometry, because the instru- 
ment + atmosphere model is wavelength-dependent, and the observables are computed 
through a difference between the data and this model. The used model is time and 
wavelength-dependent, and can be of growing complexity, depending on the wave- 
length coverage and spectral resolution of the instrument. The direct consequence is 
that this differential phase contains only a fraction of the original phas e info rmation. 

This was illustrated in lMilloud (120061 . pp. 58). lMillour et al.l(l2008l) . and lSchmitt etal 



(120091) . where a Taylor-Expansion of the observed phase <p m .(cf) as a function of wave 



number or = was presented to illustrate the information kept into the differential 
phase, and the one discarded by the process. Indeed, the observed phase, in the absence 
of nanometer-accuracy metrology, is affected by an unknown optical path difference 
(OPD) term 6jj(t), and higher order terms, which can come for example from the water- 
vapor and dry air chromatic dispersion S{j(f) : 

<f%(t, a) = f.j(<r) + 2n5ij(t)o- + Ins^o 2 + ... + o(o*) (1) 

where o{cr n ) is a negligible (or time-constant) term compared to all the other ones. 
On the other side, one can perform a Taylor-Expansion of the object phase at the obser- 
vation time: 

<p*j(cr) = a* Q + a\cr + a\o- 2 + ... + 6tp*j(cr) (2) 

where 6<p*j(cr) is the differential phase. One can see here the similarities in the 
terms in in Eq. Q] and |2] The differential data analysis process make disappear the 
l-to-« th first terms of this Taylor-Expansion, n depends on the wavelength coverage 
and spectral resolution of the instrument. The different data processing processes are 
stacking the different contributions that can be modeled, depending on the spectral 
resolution. The higher the spectral resolution, the smaller will be n, because the smaller 
the variations of o~, making more and more of the higher-order terms negligible in the 
process. 

3.1. Getting the differential phases in practice 
3.1.1. High spectral resolution 

This is the case of visible instruments (GI2T, VEGA), which have been designed and 
used with spectral resolutions of R=30,000, or the case of AMBER used in its highest 
spectral-resolution mode (R= 12,000). 

In the case where the spectral interval of interest AA is very small compared to 
the wavelength of interest A, i.e. ^4 «; 0.1, the dominant atmospheric wandering of 
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the fringes, due to the term 2n8ij{t)cr of Eq. [TJ can be considered as a constant as a 
function of cr. Therefore, that wandering is only an ensemble-movement of the fringes 
(se e Figure Q]), identical at all wavelengths. Therefore, the initial technique proposed 
bv lBeckersl(fl983h can be applied as-is, as its role is to remove the variable part of the 
phase. 
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Figure 1 . Left: Phase of the AMBER fringes as a function of time, for a bright 
object (same as in Figure 0. Right: Phase of AMBER when no fringes at all are 
present. 

Usually, the algorithms used to involve the computation of a cross-spectrum on 
each recorded frame, based on Fourier-Transforms or other algorithms (in the case of 
AMBER), and an ensemble averaging. At the very end of the processing, a cr-dependent 
residual slope of the phase may be subtracted, depending how good the white-light 
fringe was centered at observation time. 

3.1.2. Intermediate spectral resolution 

This is the case of several instruments such as AMBER, used in its medium-spectral 
resolution mode (R=1500), Keck-I in High-resolution mode (R=1000), or MIDI with 
high-spectral resolution (R=300). 

In this intermediate case (^ - 0.1), more effects become non-negligible, like 
the cr-dependence of the term 2ndjj(t)cr, which imposes a calculation of 6jj(t), and a 
frame-by-frame correction of this term (see Figure |2]). 

Therefore, for these instruments, a real-time correction (Keck, VLTI with FINITO), 
and/or a post-proce ssing correction of t his term are inclu ded in the observation/data 
processing pipeline (Woill ezet al.ll2012l : iTatulli et alj|2008h . 



3.1.3. Low spectral resolution 

This is the case of NPOI (R=38), AMBER in its low-spectral resolution mode (R=35), 
or MIDI (R=30). In this case, M > 0.1. 

In addition to the two terms described above, additional terms, like the chromatic 
dependence of the OPD, or the variation of water vapor content, cannot be neglected 
anymore. Tackling and correcting for these effects is still today an active researc h topic 
( Tubb s et al.ll200l IColavita et all2004l : IVannier et al.ll2006l : iMillour et al.ll2008l) . 
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Figure 2. "Fringe-tracks" of AMBER/VLTI, showing the fringe wandering as a 
function of time (vertical axis). The plots show the fringe power as a function of 
optical path difference (piston, horizontal axis) for a bright star, shown here for the 
three baselines. The fringe wandering range is of the order of 10's of microns over a 
few seconds, with excursions as large as 50 to 100 microns. 



One can however, for an acceptable accuracy in many cases (2-3 degrees), tackle 
the slowly- variable part of this water- v apor + dry air term by pre-determ i ned dispersio n 
laws that can be found in the literature (ICiddori[i 996: Va nnier et al.ll2006l : lMath aril2007h . 
Using the ambient conditions (temperature, pressure, humidity fraction, partial pressure 
of CO2), measured simultaneously as the observation, and the delay lines position, one 
can compute the corresponding chromatic phase effect at the time of observing. An 
example of such correction can be found in Figure [3] 



4. What does it mean? 

Differential phase is often thought to be either excessively simple, or too difficult to 
interpret. This matter of fact has several origins, which comes both from the poor 
understanding of the origin of that phase measurement (which is detailed in the previous 
section), and the fact that interpretation difficulty depends on how much th e objec t 
is resolved by th e interferometer. Ind eed, as was shown in IChelli & Petrovl (119951) : 



Lachaumel ( 2003 ) : IVannier et all (120061) . the relation between the differential phase and 



the photocenter of the object of interest is linear when the object is unresolved, or barely 
resolved (i.e. visibility is close to 1). Therefore, a direct comparison between spectro- 
interferometry and spectro-astrometry can be done in the case the object is unresolved. 

One can then get a simple qualitative interpretation of their dataset, just based 
on its morphology. One knows alr eady the "S"-shaped phases, that can be attributed 
to rotating disks dBerio et ai1ll999h . but other possibilities appeared in the meantime, 
leading to a real "zoology" of differential signal: 
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Figure 3. From top to bottom and left to right: K-band raw differential phases 
from AMBER/VLTI with intrinsic noise less than 0.05 radians, for 34 full nights 
of AMBER operation (2008). They are offset by the B x sin(z) delay function (in 
decimeters) to exhibit the dependence of the phase curvature as a function of alti- 
tude. The following graph show the curvature (2 nd order derivative) of the differ- 
ential phases as a function of delay. Th e corrected differential pha ses are shown 
afterwards, for dry air-only ICiddorl d!996lh and wet air lCiddorl (Q996) corrections. 
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the c lassical "S"-shaped phase, associated to a "W" visibility (Mei lland et al. 



2012) in a line, is linked to a rotating disk, 



on-sided phase, associated with a visibility drop in a line , can b e associated to an 
asy mmetr ic outflow or envelope lOhnaka et all (see e.g. l2009h : iBerio et all (see 



e-g- 



2oTTh . 



"double S" phase, or "W" phase, related to visibi lity variations in a line, can be 
interp reted with a bipolar outflow (see examples in lWeigelt et al.l2007l : IChesneau et al 
201 lh . 



In all other cases, where the object is resolved by the interferometer, there is no 
simple relation between the phase and the photocenter shift, though in some cases (e.g. 
Be star s or B|e1 stars), a "phase reversal" (a change of phase variation, or a phase jump, 
see e.g. Chesneau et_alj|2007|; |Meilland et al.ll2012l) can be attributed to a visibility lobe 
change. Therefore, only a proper computation of the Fourier terms calculated on a 
model, or on an image, can help to interpret the datasets. The "cookbook of differential 
phase modeling" could contain the following steps: 

• Get your freshly measured differential phase in a cool place, taking care to note 
its wavelengths and baselines, 

• Compute a pile of images of your model at as many wavelengths as your differ- 
ential phase has been measured, 

• Fourier-Transform all your images, taking care of aliasing effects, 

• Interpolate carefully this FT at the exact spatial frequency of the observation, for 
each wavelength, 

• You should get a clean phase of your model for each wavelength and each base- 
line of observation, 

• Subtract as many terms as necessary from the obtained phase, using eq. [2 

• The result is your modeled differential phase, that you can compare with your 
observed one. 



5. Conclusion 

In this proceeding, I have tried to show the tremendous advance made in the last years 
on circumstellar matter with the combination of interferometry and spectroscopy. This 
unique tool has not yet shown its huge potential in getting more physical information 
about the observed objects. 



Questions 

Q: (A. Magalhaes) Can you use the differential phase observations in order to model 
the observations, in the same way, for instance, that you do for polarization and 
position angle, where you compare with a model? 
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A: Yes, indeed, the differential phase can be computed based on the model images 
at several wavelengths, and compared to the observed one. This is actually the 
way we compare the models to the observations. 

Q: (P. Nunez) Could you please explain the use of "differential closure phase"? 

A: "Differential closure phase" & differential visibilities essentially provide the same 
information as the closure phase and squared visibilities, respectively. The main 
difference resides in the uncertainties and biases which affect these observables. 
For example, in specific cases, the "differential closure phase' can have signifi- 
cantly smaller errors than the closure phase. 
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